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Abstract 

In this article we review a theoretical framework for pair produc- 
tion from strong external electromagnetic fields. We propose a numerical 
method to solve the resulting equations of motion and present results for 
both cases of spatially homogeneous and inhomogeneous electric fields. 
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1 Introduction 

Since the prediction of e + e _ pair production from vacuum in a strong electric 
field by J. Schwinger [I], extensive theoretical research has been focused on the 
non-perturbative QED processes. While earlier, the experimental observation 
was way beyond feasible, recent advances in laser technology outperform the 
expectations. With the advent of Chirped Pulse Amplification [2] and the 
creation of femtosecond and attosecond intense laserpulses, the threshold of 
pair production is being approached by the extreme field laboratories[3j. 

At the same time, the theoretical description of the phenomena has been also 
developed. For a long time, different simple, but analitically solvable models of 
the external field were investigated but it was evident, that the full description of 
a realistic field configuration is well beyond the applicability of those models pi]. 
Recently, a Wigner function-based formalism was proposed, that can describe 
arbitrary field configurations and can describe a self-consistent description of 
the time evolution of the particle — anti-particle quantum system [HI El [7]. 

Most investigations were performed in the spatially homogeneous limit, be- 
cause the calculations were beyond the capabilities of the available computers. 
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With the advent of cheap computing devices, the calculations in multiple di- 
mensions and with spatial inhomogenities has become feasible. We present a 
numerical method that is ideal for the hardware and can efficiently solve the 
complex differential equations of pair production. 

First, we briefly review the theoretical framework; next, we discuss the pair 
yield in chirped homogeneous electric field; and finally, present first results for 
an inhomogeneous field. 

2 The Dirac-Heisenberg-Wigner formalism 

In a classical system, the one particle distribution function is used to charac- 
terize the time evolution. In the quantum case the Wigner function (W) plays 
a similar role [8j. It is defined from the following density operator [C)\ 



W(x,p,t) 



-ips 



(0\C(x,s,t)\0)d 3 s, 



(1) 



where 



C (x, s, t) = exp 



1/2 ^ ^ 

-i< I A(x + As, t)sd\ 



+ -,*),*(£- -,t) 



(2) 



Here A is the vector potential of the electromagnetic field. An equation of 
motion for the Wigner function can be derived by applying the time derivate 
to W. The derivation can be found elsewhere |7]. The equation of motion 
reads: 



D t W = -^D s [ 7 °7, W] - im [ 7 °, W] - iP { 7 °7, W) 



(3) 



where the following non-local operators are introduced: 



D t = d t + e£{x, t)V 3 - — {VsVpfS^x , t) V? +...; 

- - eh 2 ■ ■ 

Dg = W s + eB{x, t) x - — (V s Vp) 2 13(x, t) x V p ^ + 

P = V + ^(V s Vp)B(x, t)xV f +.... 



(4) 
(5) 
(6) 



with S and B denoting the electric and magnetic field respectively. 
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Figure 1: Resonant transition from adiabatic tunneling to anti-adiabatic multi 
photon pair creation (left). Zoom in around the critical frequency (u c ) for 
different chirp parameters (right). 



In the following, we will use the expansion of the Wigner function: W(x, p, t) - 
l/4[ls + «7 5 p + 7 M v M + 7 /i 75a /i + a^t^}. In this basis, the particle density can 
be calculated as: 

,(*,,) = 2 + -=±^ (7) 
ym^ + p z 

In the limit of B = and £(x,t) = £(t), one restores the well-known quan- 
tum kinetic equation j6]: 

^ = veEsJm 2 + p\/{m 2 (8) 



= ^(1 - 2f)e€yjm 2 + fj[m 2 + f) - 2u^m 2 + p 2 , (9) 



dv 

dT 

^ = 2vJm 2 +p 2 ; (10) 

where / is the particle density (u and v are auxiliary functions) and the 
kinetic and canonic momenta has the usual connection: p = (q±,q\\ — eA(t)). 
While this equation was known long ago, the total pair production-dependence 
on field parameters was not discussed for complex time-dependence of the ex- 
ternal field. 



3 Numerical Results 

To get as close to the experimental laser fields as possible in the quantum 
kinetic model we consider the following two models for the electric fields: a 
packet like £(t) = £q exp[— (t/r) 2 ] cos ((f) + ut) and a similar, but chirped field: 
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Pair distribution 



Charge distribution 




Z coordinate [A c ] Z coordinate [A c ] 

Figure 2: Pair density (left) and Charge density (right). The density scale is 
arbitrary. 

£{t) = Sq exp[— {t/r) 2 ] cos (0 + cut + ct 2 ) where the field strength is measured 
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in £ c = units, also we set h = c = 1. The system exhibits a resonant 

transition at the frequency co c = ^ between two qualitatively different regions: 
tunneling and multi photon pair creation. Figure 1 has been calculated for the 
packet-like field with Sq = 0.5 ■ S c , <fi = and r chosen in such a way as to 
have 5 cycles in the 0.1 • S max envelope. We investigated the effect of adding 
chirp to the field (Figure 1., right panel). An increase of the chirp results in the 
widening of the peaks and the rise of the average yield. If experiments aim to 
investigate this region for resonances, the chirp may become a natural control 
parameter. 

In the spatial inhomogeneous case, no simplifications are possible, and we 
must use all 16 components of the Wigner function. Problem arises because of 
the high dimensionality and the derivative expansion of the non-local operators 
(|4])-([6]). We successfully used spectral expansion to solve the system of equa- 
tions in 1+1 dimensions at the leading order of the derivative expansion. The 
method can be easily generalized to higher dimensions and can reliably resolve 
higher order terms in the derivative expansion. Details of the method will be 
discussed in a forthcoming paper. Here we illustrate the results for the following 
simple field: S(z,t) = £o ex P [ — {z 2 'jo 2 + t 2 /r 2 )]. Figure 2 show the pair density 
and the charge density of the system after the pulse. The spatial inhomogenity 
creates interesting structures in the distributions. With such methods, it be- 
comes possible to understand the microscopic spatiotemporal evolution of the 
system in a comprehensible way. 
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4 Conclusion 



We reviewed the Wigner formalism of pair production for homogeneous and 
inhomogeneous electric fields. For the homogeneous case we investigated the 
effect of chirp on the produced number of pairs. We found that the chirp 
parameter tunes the contrast of the resonance peaks. 

We developed a method based on spectral expansion to solve the Wigner 
equation in 1+1 dimensions. This method can be used to efficiently extend the 
analysis to the spatial structure of the quantum distributions characterising the 
pair production. This opens new possibilities to study more realistic laser fields 
and to provide more adequate predictions. 

This work was supported by the Hungarian OTKA Grant NK 77816. 
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